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Spin resistivity R has been shown to result mainly from the scattering of itinerant 
spins with magnetic impurities and lattice spins. R is proportional to the spin-spin 
correlation so that its behavior is very complicated near and at the magnetic phase 
transition of the lattice spins. For the time being there are many new experimental 
data on the spin resistivity going from semiconductors to superconductors. Depend- 
ing on materials, various behaviors have been observed. There is however no theory 
so far which gives a unified mechanism for spin resistivity in magnetic materials. 
Recently, we have showed Monte Carlo results for different systems. We found that 
the spin resistivity is very different from one material to another. In this paper, we 
show for the first time how the dynamic relaxation time of the lattice spins affects 
the resistivity of itinerant spins observed in Monte Carlo simulation. 

PACS numbers: 05.60.Cd Classical transport ; 75.47.-m Magnetotransport phenomena; ma- 
terials for magnetotransport ; 75.10.Hk Classical spin models ; 05.10.Ln Monte Carlo meth- 
ods 

I. INTRODUCTION 

The resistivity is an important subject in condensed-matter physics. It has been studied 
experimentally and theoretically already with old classical physics. However, only from 
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the fifties, with notions borrowed from microscopic modern physics that the resistivity has 
been viewed as a consequence of microscopic mechanisms which govern physical behaviors 
of materials through which conduction electrons travel. In this paper, we are interested in 
the resistivity caused by magnetic scattering of itinerant electronic spins by localized lattice 
spins in magnetic materials (ferromagnets and anfiferromagnets). The resulting resistivity 
is called hereafter "spin resistivity" which is to be distinguished from the resistivity due to 
spin-independent scattering, for example by phonons and non magnetic impurities. 

The spin resistivity, R, was shown to depend on the spin-spin correlation in ferromagnetic 
crystals by de Gennes and Friedel 1 , Fisher and Langer^ among others. At low temperatures 
(T), the spin- waves are shown to be responsible for the T 2 behavior of the spin resistivity in 
ferromagnets 3 ' 4 . Note however that in these calculations the itinerant electrons have been 
considered as free electrons interacting with the lattice spins, but there is no interaction 
between them. We have showed^ that when an interaction between itinerant electrons is 
introduced, the itinerant electrons can be crystallized at low T giving rise to an increase of 
R as T — > 0. Experimental data in various materials show this behavior- - — , but we would 
warn that there may be other mechanisms involved as well. At the magnetic phase transition 
temperature Tq, the spin-spin correlation diverges in magnetic materials with a second-order 
phase transition. The theory of de Gennes-Friedel predicts that R should show a divergent 
peak. However, experiments in various magnetic materials ranging from semiconductors 
to superconductors^ - — show indeed an anomaly at the transition temperature Tc, but the 
peak is more or less rounded, not as sharp as expected from the divergence of the correlation 
length. It has been shown in fact that^iIS the form of the peak depends on the length of 
the correlation included in the calculation of R: if only short-range correlations are taken 
into account, then the peak is very rounded. A justification for the use of only short-range 
correlations comes from the fact that the mean free path of itinerant spins is finite at Tc- 
When scattering is due to impurities, the peak has been shown to depend on the localization 
length^. In the case of antiferromagnets, Haas has shown the absence of a resistivity peak—. 
Our recent works using Monte Carlo (MC) simulations have shown that there is indeed an 
anomaly at Tq in various magnetic models from ferromagnets^ 1 - 1 ^, antiferromagnets^ 1 ^ to 
frustrated spin systems^ 1 ^ 4 -. The shape of the anomaly depends on many ingredients such as 
crystal structures, spin models, and interaction parameters. 

In this paper, we will show new results obtained by MC simulation when we take into 
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account the temperature dependence of the relaxation time of localized lattice spins in the 
simulation. We will show that this temperature dependence affects the shape of the peak in 
the phase transition region. 

Section [XT] is devoted to a description of the general model and the MC method. We 
introduce in this section the temperature dependence of the relaxation time. Results are 
shown and discussed in section II III for both ferro- and antiferromagnets in terms of critical 
slowing-down. Concluding remarks are given in section IIVI 

II. MODEL AND METHOD 

The model we use in our MC simulation is very general. The itinerant spins move in a 
crystal whose lattice sites are occupied by localized spins. The itinerant spins are assumed 
to be of Ising type, but the method of simulation can be used for other spin models 23 . The 
localized spins may be of Ising, XY or Heisenberg models. Their interaction is usually limited 
to nearest neighbors (NN) but this assumption is not necessary. It can be ferromagnetic or 
antiferromagnetic . 

A. Interactions 

We consider a thin film of a given lattice structure where each lattice site is occupied 
by a spin. The interaction between the lattice spins is limited to NN with the following 
Hamiltonian : 

Hi = - zJ Ji,jSi- S i (!) 

(ij) 

— * 

where Si is an Ising spin whose values are ±1, Jij the exchange integral between the NN 
spin pair Si and Sj. Hereafter we take Jij = J for all NN spin pairs, for simplicity. As a 
convention, ferromagnetic (antiferromagnetic) interaction has positive (negative) sign. The 
system size is N x x N y x N z where Ni(i = x,y,z) is the number of lattice cells in the i 
direction. Periodic boundary conditions (PBC) are used in the x and y directions while the 
surfaces perpendicular to the z axis are free. The film thickness is N z . 

We define the interaction between the itinerant spins and the localized lattice spins as 
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follows 

H r = -^IijCTi.Sj (2) 

i,j 

where Oi is the Ising spin of the i — th itinerant electron and I^j denotes the interaction 
that depends on the distance between an electron i and the spin Sj at the lattice site j. For 
simplicity, we use the following interaction expression 

h :3 = Joe-"* (3) 

where = \ fi — fj\, Iq and a are constants. In the same way, interaction between itinerant 
electrons is defined by 

Urn = X K '-J a ''- 'i ( 4 ) 

K id = Koe-fo (5) 

with Kij being the interaction that depends on the distance between electrons i and j. The 
choice of the constants K and (5 will be discussed below. 

Dynamics of itinerant electrons is ensured by an electric field applied along the x axis. 
Electrons travel in the x direction, leave the system at the end. The PBC on the xy planes 
ensure that the electrons who leave the system at one end are to be reinserted at the other 
end. For the z direction, we use the mirror reflection at the two surfaces. These boundary 
conditions are used in order to conserve the average density of itinerant electrons. One has 

U E = -eel (6) 

— * 

where e is the charge of electron, e the applied electric field and £ the displacement vector 
of an electron. 

Since the interaction between itinerant electron spins is attractive, we need to add a 
chemical potential in order to avoid a possible agglomeration of electrons into some points 
in the crystal and to ensure a homogeneous spatial distribution of electrons during the 
simulation. The chemical potential term is given by 

H c = DV r n(r) (7) 

where n(r) is the concentration of itinerant spins in the sphere of D 2 radius, centered at r. 
D is a constant parameter appropriately chosen. 
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B. Simulation Method 

The procedure of our simulation can be split into two steps. The first step consists 
in equilibrating the lattice at a given temperature T without itinerant electrons. When 
equilibrium is reached, in the second step, we randomly add iVo polarized itinerant spins 
into the lattice. Each itinerant electron interacts with lattice spins in a sphere of radius D\ 
centered at its position, and with other itinerant electrons in a sphere of radius D 2 . 

The procedure of spin dynamics is described as follows. After injecting A itinerant 
electrons in the equilibrated lattice, we equilibrate the itinerant spins using the following 
updating. We calculate the energy E Q u of an itinerant electron taking into account all 
interactions described above. Then we perform a trial move of length i taken in an arbitrary 
direction with random modulus in the interval \R\^R 2 \ where K\ — and R 2 = a (NN 
distance), a being the lattice constant. Note that the move is rejected if the electron falls in 
a sphere of radius r centered at a lattice spin or at another itinerant electron. This excluded 
space emulates the Pauli exclusion. We calculate the new energy E new and use the Metropolis 
algorithm to accept or reject the electron displacement. We choose another itinerant electron 
and begin again this procedure. When all itinerant electrons are considered, we say that 
we have made a MC sweeping, or one MC step/spin. We have to repeat a large number of 
MC steps/spin to reach a stationary transport regime. We then perform the averaging to 
determine physical properties such as magnetic resistivity, electron velocity, energy etc. as 
functions of temperature. 

We emphasize here that in order to have sufficient statistical averages on microscopic 
states of both the lattice spins and the itinerant spins, we use the following procedure: after 
averaging the resistivity over Aq steps for "each" lattice spin configuration, we thermalize 
again the lattice with A 2 steps in order to take another disconnected lattice configuration. 
Then we take back the averaging of the resistivity for Aq steps for the new lattice config- 
uration. . We repeat this cycle for A3 times, usually several hundreds of thousands times. 
The total MC steps for averaging is about 4 x 10 5 steps per spin in our simulations. This 
procedure reduces strongly thermal fluctuations observed in our previous work 22 . 

Of course, the larger A x and A 3 are the better the statistics becomes. The question is 
what is the correct value of Ni for averaging with each lattice spin configuration at a given 
T? This question is important because this is related to the relaxation time -r^of the lattice 
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spins compared to that of the itinerant spins, tj. The two extreme cases are i) t l ~ tj, one 
should take Ni — 1, namely the lattice spin configuration should change with each move 
of itinerant spins ii) tl ^> ti, in this case, itinerant spins can travel in the same lattice 
configuration for many times during the averaging. 

In order to choose a right value of N±, we consider the following temperature dependence 
of tl in non frustrated spin systems. The relaxation time is expressed in this case as 25 - 

- - 

Tl ~ |1 - TjT c \ zv 

where A is a constant, v the correlation critical exponent, and z the dynamic exponent. From 
this expression, we see that as T tends to Tc, tl diverges. In the critical region around Tc the 
system encounters thus the so-called " critical slowing down" : the spin relaxation is extremely 
long due to the divergence of the spin-spin correlation. In our previous papers^&^i^ 1 ^, we 
did not take into account the temperature dependence of T£. We propose to study here the 
spin resistivity using Eq. flH]). 
We define spin resistivity p as : 

p = - (9) 

n e 

where n e is the number of itinerant electron spins crossing a unit slice perpendicular to the 
x direction per unit of time. 



C. Choice of parameters and units 

The spin resistivity is dominated by the two interactions Eqs. ([2]) and As said 
earlier, our model is very general. Several kinds of materials such as metals, semiconductors, 
insulating magnetic materials etc. can be studied with our model, provided an appropriate 
choice of the parameters. For example, non magnetic metals correspond to ijj = Kij = 
(free conduction electrons). The case of magnetic semiconductors corresponds to the choice 
of parameters K and Iq so as the energy of an itinerant electron due to the interaction 
l-l r should be much lower than that due to % m , namely itinerant electrons are more tightly 
bound to localized atoms. Note that ri m depends on the concentration of itinerant spins: for 
example the dilute case yields a small T-L m . We will show below results obtained for typical 
values of parameters which correspond more or less to semiconductors. The choice of the 
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parameters has been made after numerous test runs. We describe the principal requirements 
which guide the choice: i) We choose the interaction between lattice spins as unity, i. e. 
| J\ — 1, ii) We choose interaction between an itinerant and its surrounding lattice spins so 
as its energy Ei in the low T region is the same order of magnitude with that between lattice 
spins. To simplify, we take a = 1. This case corresponds more or less to a semiconductor, 
as said earlier, iii) Interaction between itinerant spins is chosen so that this contribution 
to the itinerant spin energy is smaller than Ei in order to highlight the effect of the lattice 
ordering on the spin current. To simplify, we take /3 = 1, iv) The choice of D is made in 
such a way to avoid the formation of clusters of itinerant spins (agglomeration) due to their 
attractive interaction [Eq. (J^D], v) The electric field is chosen not so strong in order to avoid 
its dominant effect that would mask the effects of thermal fluctuations and of the magnetic 
ordering, vi) The density of the itinerant spins is chosen in a way that the contribution 
of interactions between themselves is much weaker than E{, as said above in the case of 
semiconductors. 

A variation of each parameter respecting the above requirements does not change qual- 
itatively the results shown below. Only the variation of D\ in some antiferromagnets does 
change the results (see Ref.-). 

The energy is measured in the unit of |J|. The temperature is expressed in the unit of 
| J\/kB- The distance (Di and D 2 ) is in the unit of a. 

III. RESULTS AND DISCUSSION 

In this section, we show for comparison the results obtained with and without temperature 
dependence of the lattice relaxation time for both ferromagnets and antiferromagnets. In 
each case we use the same set of interaction parameters in order to outline the effect of the 
temperature-dependent relaxation time. 

In this paper we use the lattice size N x = N y = 20 and N z = 8 and we consider the body- 
centered cubic (BCC) lattice for illustration. The lattice constant is a. The spin resistivity 
is calculated with N = (N x x N y x N z )/2 itinerant spins (one electron per two lattice cells). 
Except otherwise stated, we choose interactions Iq = 2, Kq = 0.5, D\ = a, D 2 = a, D = 0.5, 
e = 1, Nq = 1600, and ro = 0.05a. A discussion on the effect of a variation of each of these 
parameters is given above. 
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Note that, due to the form of the interaction given by Eq. (jSJ), the itinerant spins have 
a tendency to form compact clusters to gain energy. This tendency is neutralized by the 
concentration gradient term, i. e. a chemical potential, given by Eq. ((?]). The value of D has 
to be chosen so as to avoid an agglomeration of itinerant spins. This choice depends of course 
on the values of D\ and D2. Examples have been shown elsewhere.- 1 ^ For the temperature 
dependence of the lattice relaxation time tl, we take v = 0.638 (3D Ising universality) and 
z = 2.02.— By choosing A — 1, we fix tl — 1 at T — 2Tq deep inside the paramagnetic 
phase far above Tc- This value is what we expect for thermal fluctuations in the disordered 
phase. 

Figure [1] shows the spin resistivity R in a BCC ferromagnet. Note that the transition 
temperature for this thin film of size 20 x 20 x 8 with Ising spins interacting via the NN 
coupling is Tq — 6.35. Several remarks are in order: 

i) The results obtained with and without temperature-dependent relaxation time for 
T < Tc coincide with each other 

ii) At Tc, for the set of parameters used here, the results using the temperature- 
independent relaxation shows a broad maximum above Tc while those using the 
temperature-dependent relaxation strongly decreases at Tc giving rise to a sharp peak. 

We show now in Fig. [2] the spin resistivity R in a BCC antiferromagnet. Note that the 
transition temperature is the same as that of the ferromagnet counterpart shown above. 
Here we observe that R in the case of temperature-dependent relaxation is lower than that 
in the case of temperature-independent one in the whole temperature range. Note that the 
value of the peak is much smaller here than in the ferromagnet case. 

We show in Fig. [3] the two curves of ferromagnet and antiferromagnet with T-dependent 
relaxation time. We observe here that below Tc, the resistivity of antiferromagnet is higher 
than that of ferromagnet, while for T > Tc the reverse is true. 

It is interesting to calculate the relaxation time 77 of the itinerant spins. We define 77 in 
the simulations as the MC time (in unit of one MC step/spin) between two "MC collisions", 
namely the lapse of time between two "rejections" of a spin to advance. Of course this 
quantity is averaged over all itinerant spins and over the simulation time. Figure H] shows 
rf 1 obtained by simulation using 77. As for R seen above, the temperature dependence and 
independence are markedly different only for T > Tc- The same thing is observed for the 
case of antiferromagnet shown in Fig. [5] 
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FIG. 1: BCC ferromagnetic thin film: Resistivity R with temperature-independent relaxation 
(white circles) and temperature-dependent relaxation (black circles) in arbitrary unit versus tem- 
perature T, in zero magnetic field, with electric field e = 1, Jq = 2, Kq = 0.5. 




FIG. 2: BCC antiferromagnetic thin film: Resistivity R with temperature-independent relax- 
ation (white circles) and temperature-dependent relaxation (black circles) in arbitrary unit versus 
temperature T, in zero magnetic field, with electric field e = 1, Iq = 2, Kq = 0.5. 

Figure [6] shows rf 1 for both ferro- and antiferromagnetic cases, for comparison. The 
antiferromagnet has rf 1 larger at T < To- Note that the resistivity is proportional to rf 1 . 

Let us discuss the reason why the temperature dependence of the lattice relaxation time 
affects so strongly the shape of the spin resistivity for T > Tc- First, we emphasize on two 
"empirical" rules that we observed and verified in a number of cases: 

• a) itinerant spins move easily when they are energetically unstable. The electric field 
then drives them easily forward. On the other hand, when they are "at ease" with 
surrounding spins, namely their energy is low, they will not move easily. We have 
checked this rule by calculating their velocity as a function of their energy^ 
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FIG. 3: BCC ferromagnetic and antiferromagnetic films: Resistivity R with temperature- 
dependent relaxation for ferro- (black circles) and antiferromagnet (white circles) in arbitrary 
unit versus temperature T, in zero magnetic field, with electric field e = 1, Jo = 2, Kq = 0.5. 
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FIG. 4: BCC ferromagnetic film: Inverse of relaxation time of itinerant spins tJ v calculated with 
temperature-independent (white circles) and temperature-dependent (black circles) of the lattice 
spins versus temperature T, in zero magnetic field, with electric field e = 1, Iq = 2, Kq = 0.5. 

• b) in the case where the energy of an itinerant is low, the move of itinerant spins 
depends on the energy difference between its initial and final positions. Consider 
the ordered phase of the lattice: the energy at any point is very low and the energy 
difference between any two points is close to zero (ordered state). So, by the MC 
updating criterion, the electric field dominates again the move of itinerant spins. This 
explains the very small resistivity at low T with respect to that at high T (except 
when T — > where other mechanisms come to play). 

For the effect of Tl, several important points are in order: 

i) For T < T c the lattice is ordered, therefore itinerant spins do not see the difference 
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FIG. 5: BCC antiferromagnetic film: Inverse of relaxation time of itinerant spins rf 1 calculated 
with temperature-independent (white circles) and temperature-dependent (black circles) relaxation 
time of the lattice spins versus temperature T, in zero magnetic field, with electric field e = 1, 
I = 2,K = 0.5. 




FIG. 6: BCC ferromagnet and antiferromagnet cases: Inverse of relaxation time of itinerant spins 
rf 1 calculated with temperature-dependent relaxation time of the lattice spins for ferromagnet 
(black circles) and antiferromagnet (white circles) versus temperature T, in zero magnetic field, 
with electric field e = 1, Iq = 2, Kq = 0.5. 

when the lattice changes its microstates more frequently or less frequently. This explains the 
same values obtained for R with and without temperature dependence of tl in ferromagnets. 
In antiferromagnets, one observes a small difference due to the presence of lattice down spins 
which act differently on up-polarized itinerant spins. 

ii) For T > T C) the lattice is disordered: the lattice spins are frequently nipped. Itinerant 
spins have to move constantly to accommodate themselves to the fluctuating environment. 
Thus, Tj is long by definition because there are very few rejections to move. Consequently, 
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R is small in the paramagnetic phase. 

iii) Finally, it is striking to observe a strong correlation between 77 and 77: Since 77 is 
very large in the transition region where the lattice is in the regime of critical slowing down, 
itinerant spins have time to find themselves in energetically favorable positions. Once they 
are there they refuse to move (first rule mentioned above). As a consequence 77 is very small 
(for example 77 = 1 if they refuse to move at every update trial). R is thus very high. We 
have showed the inverse of 77 in Figs. HI and [H] because R is inversely proportional to 77. 
The correlation between 77 and 77 is thus "high 77 corresponds to low 77" and vice- versa. 

Now, let us show the effect of the choice of A of Eq. (JSJ) in Fig. [7J The higher A (i. 
e. higher 77) induces an increase of R near Tc but gives the same value as T is far away 
from the critical point. Thus, the width of R at the transition temperature can serve as a 
measure of the relaxation time of the lattice spins. 
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FIG. 7: BCC ferromagnet. Effect of the choice of the constant A: i) A = 1 (black circles), ii) 
A = 2 (white circles) versus temperature T, in zero magnetic field, with electric field e = 1, Iq = 2, 
Kq = 0.5. See text for comments. 



IV. CONCLUSION 

In this paper, we have shown the effect of the temperature dependence of the relaxation 
time on the spin resistivity for both ferromagnetic and antiferromagnetic films with BCC 
lattice structure. 

In the ferromagnetic case, the long relaxation time in the critical region compared to 
that of the paramagnetic phase gives rise to a sharp peak of the spin resistivity at T c . 
The resistivity in the low-T region is insensitive to the relaxation time while in high-T 
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region, the resistivity is much smaller than that obtained with the temperature-independent 
relaxation time. The same tendency is observed for the antiferromagnetic case: while the 
spin resistivity in the case of temperature-independent relaxation time does not show a 
peak at Tc, the extremely long relaxation in the critical region with respect to that of the 
paramagnetic phase gives rise to a pronounced rounded peak at Tc- It is very interesting to 
study other systems such as spin glasses where the relaxation time is extremely long even 
at temperatures far below T c . 
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